THE VARIANCE OF THE SHOCK IN THE HAD PROCESS 



CRISTIAN F. COLETTI, PABLO A. FERRARI AND LEANDRO P. R. PIMENTEL 



Abstract. We consider the Hammcrsley-Aldous-Diaconis (HAD) process with sinks and 
sources such that there is a microscopic shock at every time t; denote Z{t) its position. We 
show that the mean and variance of Z{t) are hnear functions of t and compute cxpUcitely 
the respective constants in function of the left and right densities. Furthermore, we 
describe the dependence of Z{t) on the initial configuration in the scale and, as a 
corollary, prove a central limit theorem. 



1. Introduction 

Let S and W be one-dimensional Poisson processes and let P be a two dimensional 
Poisson process of rate 1. Assume they are homogeneous and mutually independent. The 
Hammersley-Aldous-Diaconis process [6l [1], or shortly the HAD process, T-C{S,W, P) = 
{Hs, s G [0,t]) has been constructed by Groeneboom [2] in the square [0,x] x [0,t] as a 
deterministic function of S, W and P as follows. The point configuration Hg represents 
the position of particles. At time zero the particles start at the points in S, called the 
sources. Then (see Figured]), at the first s > such that {y,s) is in P for some y G [0,x] 
or s is in W, the closest particle to the right of y jumps to y if {y, s) is in P or to if s is 
in W (points in W are called sinks). If there is no particle to the right of y, then a new 
particle is added at y at time s. The new configuration does not move until the second s 
such that {y, s) is in P for some y or s is in W, when the second jump occurs, and so on 
until time t. In other words, we define the process inductively as follows: Hq = S and for 
s > 0, 



Hs^ ii s VJ {s', (y, s) e P for some y' G [0, x]} 

[h,. \ {R{y, Hs-)}] U {y} if (y, G P U {(0, s'), s' G W} 



where R{y,H) = inf{y' G i/ U {x} : y' > y} for H C (0,x). Let N be the positions of 
particles at time t, and let E be the set of times a new particle enters the system from the 
right through the vertical axis {x} x [0,t). 

Let A > and p > 0. Assume that S has intensity A, W has intensity p and P has 
intensity 1. Define the point process S' by removing the first sink of S and consider the 
corresponding HAD process Hi^S', W, P). For time t > the coupled processes 7^(5", W, P) 
and 7i(S", VF, P) will differ at one point denoted Z{t) and called a second class particle 
(Figure [1]). 
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S (".0) 



Figure 1. The black points represent the Poisson points S (sources), W 
(sinks) and P. The polygonal lines represent the trajectory of HAD particles 
while the trajectory of the second class particle is represented by the dashed 
polygonal line. 



Results. Our main results are the computation of the mean and variance of Z{t) and the 
dependence of fluctuations on the initial configuration. 



Theorem 1. Assume that Xp > 1. Then for all t > we have 



E{Z{t)) 



P 
A' 



t. 



Let 



Var(Z(t)) = 4^ 



RiS)^l{p-\) andi?(iy)^i(p-l). 



;i.2) 



In the next theorem we show that in the scale ^/t the position of Z{t) is the same as the 
following random variable, depending only on 5* and W: 



I [0, tR{S)] ns\ + \ [0, tR{w)] n w\ 

where \A\ denotes the number of points in the (finite) set A. 
Theorem 2. Assume that Xp > 1. Then 



E 



lim 

f— »oo 



{zit)-Nit)y 



0. 
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and 



t— >oo 



lim 



Z(t)-E(Z(t)) 
v/Var(Z(t)) 



A/'(0, 1) in distribution, 



where A/'(0, 1) a standard normal random variable. 

Brief historical overview. The HAD process has a genealogical relation with the famous 
Ulam's problem, which studies the limit behavior of the longest increasing subsequence L„ 
in a random permutation of the numbers 1, . . . , n. At the beginning of the seventies Ham- 
mersley [6] proposed to solve this problem as follows. Assume . . . , are n 

independent random points distributed uniformly in the rectangle [0, x] x [0,t]. These 
points specify a uniform random permutation vr by setting the point with the z'th smallest 
t-coordinate has the 7r(z) smallest x-coordinate. The length of the longest increasing sub- 
sequence of TT equals the maximal number M{x, t) of points on an up-right path from (0, 0) 
to (x, t) (last-passage percolation). The variable M{x,t) has the so called super-additivity 
property [7], which was the key to show that \fnLn — > c. Hammersley conjectured that 
c = 2 by presenting a non-rigorous but quite reasonable argument based on hydrodynam- 
ical local equilibrium for a related discrete-time particle system. Later Logan and Sheep 
[To] and Vershik and Kerov [11] proved c = 2 with combinatorial methods. 

In the middle of the nineties, Aldous and Diaconis [1] constructed a continuous time 
version of the Hammersley process and use Rost [12] coupling ideas, developed in the 
exclusion process context, to show that c = 2. This is the reason we call this process 
Hammersley-Aldous-Diaconis process. Groeneboom [9] introduced the construction in the 
quadrant we use here. In this setup M(x,t) equals the number of particles in [0,x] x \t\ 
plus the number of particles in {0} x [0,t] (we are allowing an up-right path to pick 
Poissonian points in the horizontal or vertical lines). Thus, the trajectory of the n'th 
particle corresponds to the boundary of the region {M(x, t) < n — 1} (Figure 1). 

The macroscopic hydrodynamical behavior of the HAD process is given by the Burgers 
equation dtu -|- d^giu) = 0, with g{u) = ^ [13]. When the initial conditions are u{r, 0) = A 
for r > and u(r, 0) = ^ when r < 0, the solution has a drastic changing when the product 
Xp goes from values smaller than 1 to values greater than 1. For Ap = 1, the system is 
stationary and the solutions are constants in time. For Ap < 1 it develops a rarefaction 
front while for Ap > 1 it develops a (macroscopic) shock. Since our results concern the 
shock regime, from now on we assume Ap > 1. 

The microscopic structure of a shock is described by the so called second-class particle. 
Our main result is the explicit expression (11. 2p for the variance of the second class particle 
at any given time t. The method, inspired in work by one of the authors and Fontes [H E] 
for the exclusion process, relates the mean and the variance of the flux of particles through 
the origin with the mean and the variance of the position of a single second class particle by 
a linear equation. In this way, if we calculate one, we get the other. However those works 
show that for the exclusion process the variance of a second class particle is asymptotically 
linear in time; the constant can also be computed. In the HAD process context, we show 
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how to refine this method to get an exact formula for the variance of a second class particle 
at any given time t. 

2. Preliminaries, coupling and Burke's theorem 
Notation. For a generic point process M we denote 

^ f{y,M)Midy):= J] /(x,M) (2.3) 

y£Mr\[0,x] 

and recall the Slivnyak-Mecke Theorem: if M is a Poisson process with intensity A, then 

f{y,M)M{dy))= Ef{y, M)X{y)dy (2.4) 

The stationary regime and the Burke's theorem. Let 7 > and and be 

Poisson processes with intensity 7 and 7"^ respectively. Then x = '^{^x^^xy ^) ^ 
stationary HAD process with parameter 7. Cator and Groeneboom [2] proved the following 
Burke's theorem: A'^ and are independent Poisson processes with the same density as 
and W^, respectively. 

Coupling two stationary process. Let Xp > 1. Let S*^, W^r, I and J be mutually 
independent one-dimensional Poisson point processes with intensities ^, j, X — ^ and p — j 
respectively. Set 

S^ = Sr, + I andWr, = W^ + J. (2.5) 

Thus, 5*0- and are two (independent) one-dimensional Poisson processes with intensity 
A and p respectively. Let P be the rate-1 two-dimensional Poisson process mentioned in 
the Introduction. 

We run three HAD processes simultaneously with different boundary conditions but with 
the same P. This is a realization (</>, a, rj) of the so called basic coupling defined by 

<P = niS^, Wr,, P),a = niS^, W,, P) and t] = n{s^, W^, P) . 

For X ^ ^6 denote by xi^it) the number of x-particles (particles that count for 

the process x) the interval [0,x] at time t, with the convention that particles escaping 
through a sink during the time interval [0, t] are located at 0. The flux of the discrepancies 
between a and 77 through the space-time line {0,0)-{x,t) is the process defined by 

^{x,t) = a{x,t)-7]{x,t). (2.6) 

There are two kinds of discrepancies: i) those starting at a point y at time in / (in 
the X-axis); ii) those starting at a added sink in a time s in J (in the t-axis). The position 
at time t of the second class particle starting at time zero at site y is denoted Z(y^o)(^); or 
shortly Zy{t), and the position at time t of the second class particle starting at site at 
time s is denoted Z(o,s)(t), or shortly Zs(t). 
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The flux of second class particles defined in (12.61) is also given by 

e(x,t) = #{ye/ : Z(,,o)(t) <x}-#{se J : Z(o,.)(t) >x}. 

:= Ux,t)-^_{x,t) (2.7) 

3. The mean and the variance of the flux of second class particles 

In this section we show the following. 
Proposition 1. For all x,t >0 

E(e(x,t)) = (A-l)a;- (p-^)t. (3.8) 

Var(e(a:, t)) = (^A - + (p - . (3.9) 

Denote by M := |M|, the total number of points of a generic finite point process M and 
X = xi^^t)- With this notation, 

X = N^ + W^ = S^ + E^, (3.10) 

for X £ v} ^i^d 

e = ^-rj. (3.11) 
Proof of Proposition [1], By definition (12.61) and Burke's theorem, 

E(e(x,t)) = E(^) = E{a) - E(rj) = E(iV, + W^- {N, + W,)) = (Ax + A) - (- + pt) . 

A p 

This shows (13.81) . 

To compute the variance combine first (12. 5p and Burke's theorem to get 

Cov{S^,Wrj) = Cov{E^,Nr,) =0. (3.12) 

Hence, by f lXTUD . 

Cov(a,r/) = CoY{S„ + E^,Nr, + Wr,) 

= Gov [Sr, + /, A^,,) + Gov [E^, + J) 

= Cov{S^,Nr,)+Cov{E,,W^). (3.13) 
Analogous (and simpler) reasoning for x £ {v^ gives 

Var(x) = Gov(E^, W^) + Gov(iV^, S^) . (3.14) 

On the other hand, since x = + 

Var(x) = Var(iVx) + Var(iy^) + 2Gov(iVx, W^) (3.15) 
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Using = Sy_ + E^ — N-^ and independent of A^^, 

Cov(iV,, W^) = Cov(iV„ S^ + E^- N^) = Cov(iV„ S^) - Var(iV,) . 
which imphes 

Var(x) = Var(iy^) - Var(iV^) + 2Cov(iV^, S^) . (3.16) 
Identities (131H and (IHTTHD imply 

Cov{E^, W^) - Cov{N^, 4) = Var(iy^) - Var(iV^) (3.17) 

Finally we compute Var(^) using (I3.13P and (13.141) : 

Var (^) = Var (a) + Var {f}) - 2Cov {a, if) 

= Cov{E^, Wr,) - Cov{N,„ S^) + Cov(iV,, S^) - Cov{E^, W^) 

Var(iy^) -Var(A^^)) - (Var(iy<,) - Var(A^,) 



= {pt-^)-{j-Xx), (3.18) 
where in the third line we used (I3.17p . □ 

4. The mean and the variance of a second class particle 

Lemma 1. For al x,t > 0, 

E^{x,t)+ = {X--)\x- [ ¥{Z{t) > z)dz] . (4.19) 
P ^ Jo ^ 

E^{x,t)- = {p-j)^^t-j^F{Z{u)<x)du^ (4.20) 

Proof of Lemma [1], By definition (12.71) . 

e(x,t)+ = / l{Zy{t)<x}I{dy) 
Jo 

where Zy(t) denotes the position at time t of a second-class particle that starts from the 
discrepancy at ?/ G [0,x]. By translation invariance and (12. 4p . 

Ee(x,t)+ = (A--) / ¥{Z{t)<x-y)dy. (4.21) 
P Jo 

= (X--) [ F{Z{t) < z)dz 
P Jo 

which proves (14.191) . Analogously, (14.201) follows from 

E^{x,t)- = E(^j F{Zs{t) > x)J{ds)y (4.22) 

□ 
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The main step to calculate the mean value of Z is the following. 
Proposition 2. For all x,t >0 



[ F{Z{t) > z)dz = ^ [ F{Z{u) <x)du. (4.23) 
Jo Jo 



Proof of Proposition [H Combining fICTD . fOOj) with fM . 

(A-i)a;- (p-i)t = E(e(x,t)) = E(e(a;,t)+) -E(e(x,t)_) 

= (A-i){a;-£p(Z(t)>^)rf.} - (p-i)|t-£p(Z( 



u) < X 



that is, 



(A--) / P(Z(t)>2)rf^ = (p--) / P(Z(M)<x)rfM (4.24) 
P Jo A Jo 

which is equivalent to fl4.23p for Ap > 1. □ 

Proof of (11.11) . It follows directly from Proposition [21 

E(Z(t)) = lim / P(Z(t) > = lim ^ / F(Z{u) <x)du = ^t. 

x^oo Jq x^oci X Jq a 

□ 

A similar idea works to compute Var(Z'(t)): 
Proposition 3. For all x,t > 

{p-j)(t + j^F{Z{u)<x)du^ (4.25) 
= (A-^)^(^2^ zF{Z{t)> z)dz- (^j^ P(Z(t) > + Var(e(x,t)_) 

Proof of (11.21) . Taking x — oo in both members of ( ]4.25p . 

{p-\)2t = (A-i)Var(Z(t)), 

from where (11.21) follows. (Notice that, for fixed t > 0, ^-(x, t) decreases to when x oo.) 

□ 

Proof of Proposition [S]. Noticing that ii y < y and Zy{t) < x then Zy{t) < x, 

l{Zy{t) < x}l{y ^ y}I{dy)I{dy) (4.26) 
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as we do for the square of a sum of functions with values in {0, 1}. The expected value of 
the double integral in the above equation is given by 

l{Zy{t)<x}l{yy^y}I{dy)I{dy)) 

= 2{\--y r !\{Z{t)<x-y)dydy 
P Jo Jo 

= 2{\--Y [ yF{Z{t)<x-y)dy, 
P Jo 

= 2{\--^^[x F{Z{t) < z)dz- z¥{Z{t) < z)dz^ 

= 2[\--^^[x^-xj^ F{Z{t)> z)dz-^ + z¥{Z{t) > z)dz'^ (4.27) 
using translation invariance. Take the expectation of (14.261) using (I4.19p and (14.271) to get 
E{^{x,t)l) = {\--^{x-j^¥{Z{t)>z)dz^ (4.28) 

+ {\--Y(x^ ~2x f ¥{Z{t) > z)dz + 2 [ zF{Z{t)>z)dz 
P ^ Jo Jo 



2 

iz' 



P 

On the other hand, taking the square of ( 14.19p . 

(E^(a;,t)+)^ = (A- ^)^|x2-2x^ ¥{Z{t) > z)dz + P(Z(t) > (4.29) 

Substracting flCTj) from K28\f : 

Var(e(x,t)+) = {X--)(x- f F{Z{t) > z)dz) 

P ^ Jo ' 

+ (A-i)^(^2^ zF{Z{t)> z)dz- (^J^ F{Z{t)>z)dz 
that is, using ( 13. 9p and Var(^(x, t)) = Var(^(x, t)+) + Var(^(x, t). 

(p-i)t = -{x-^)[jy{z{t)>z)dz) 

+ (A - -)Y2 / zF{Z{t) > z)dz -( ¥{Z{t) > z)dz 
P ^ Jo ^ Jo ' 

+ Var(e(x,t). 

Use fOSll to get 

ip~j){^ + j^ P(Z(n) < x)du^ 

= (A-^)^(^2^ zF{Z{t) > z)dz- (^J^ ¥{Z{t) > z)dzy^ + VaT{^{x,t). 
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which shows fl4.25p . □ 



5. The dependence on the initial condition 

Recall that 

1 l-tRiS) ptR{W) 

N{t) = (A - -) Cot - ( / S{dz) + / W{ds) 

here 

Co = 3(p-^), R{S) = Up-^^ andi?(iy) = -(p-^ 



Thus, 

E{N{t)) = E{Z{t)) = and Var(A^(t)) = Var(Z(t)) = Dt (5.30) 

A 

where 

2(p-i) 



D 



(by calculating N{t) and using Theorem [1]) . The next step is to prove: 
Lemma 2. 

^.^^ Cov(Z(t),iV(t)) ^ 

Proof of Lemma [2l The proof follows in the same lines of an analogous result proved by 
Ferrari [5] in the TASEP context. To sketch the idea, denote by Z^^^{t) (respectively, Z^^^(t)) 
the position at time t of a second class particle with respect to an initial configuration con- 
ditioned to have a particle at position z (respectively, conditioned not to have a particle 
at position z). Then we can get (by coupling) 

Cov(Z(t), / Sidz)) = -\ E{Z'^%t) - Z'^\t))dz (5.31) 

Jo Jo 

and 

Cov(Z(t), / W{ds)) = -p E{Z'^%t) - Z''^\t))ds. (5.32) 

Jo Jo 

Now, we claim that for all e > 

lim sup |E(Z"l°(t) -Z"l^(t)) - , \ I = 0. (5.33) 

*-°°.G[0,t(i?(.)-e)] (^-p) 

To sketch the idea, call (see Section [3]) Zn^{t) the position at time t of the second class 
particles starting at ^ conditioned to have one particle at z. Let Yt be the position at time 
t where the processes with initial configuration being different only at z differ by time t. 
Let Ti be the first time the trajectories of Yt and Z^^^{t) meet each other. Note that Ti 
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is finite with probability one. Before Ti,Z'\^{t) = After ri,Z^I^(t) = ^i'^(t) and 

= Z2^{t). Since 

E(zf(t)-Zf(t)) = — 1^, 
we get fl5.33l) . Thus, Lemma [2] follows by combining (15.31 1) . fl5.32p together with fl5.33l) . □ 
Proof of Theorem [H By flOOj) . 



E 



{Z{t)-N{t)y =2Dt -2CoY {Z{t),N{t)) 



and hence, Theorem [2] follows from Lemma [21 □ 
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